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Abstract 

A study of the J/ipmr and J/tpKK systems, treating them as cou- 
pled channels, has been made by solving the Faddeev equations, with 
the purpose of investigating the possibility of generation of the J^'^ = 
1 , y(4260) resonance due to the interaction between these three 
mesons. In order to do this, we start by solving the Bethe-Salpeter 
equation for the two pseudoscalar and for the vector-pseudocalar me- 
son systems using the amplitudes obtained from the lowest order chiral 
Lagrangians as potentials. With the t-matrices generated from these 
potentials, which contain the poles of the a, /o(980) and ao(980) res- 
onances for the pseudoscalar-pseudoscalar system and the pole of the 
A (3872), alongwith other new charmed resonant states, for the vector- 
pseudoscalar system, we solve the Faddeev equations. As a result, we 
get a peak around 4150 MeV with a width ~ 90 MeV when the in- 
variant mass of the two pseudoscalars is close to that of the /o(980). 

1 Introduction 

An enhancement in the data for the -k^-k^ J /ip invariant mass spectrum was 
found near 4.26 GeV by the BABAR collaboration in a study of the e^e~ — )■ 
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liSR'^^'^^ J li' process[lJ. A fit to this data set was made by assuming a 
resonance with 4.26 GeV of mass and 50 to 90 MeV of width pQ. The 
resonance was named as the y(4260) and it was found to be of J^^ = 1 
nature. Later on, an accumulation of events with similar characteristics in 
the 7T~^7i~J/ip, T^^T^^J/ip and the K~^K~J/ip mass spectra was reported by 
the CLEO collaboration [21 [3], thus confirming the results from BABAR. 
Following these works, the BELLE collaboration obtained the cross sections 
for the e~^e~ — )■ ti^tt~ J /ip reaction in the 3.8 to 5.5 GeV region [1], by keeping 
all the interactions in the final state in S-wave, and found a peak at 4.26 GeV 
and a bump around 4.05 GeV. 

Although the y(4260) does not seem to fit in to the charmonium spec- 
trum of the particle data group [5] known up to ~ 4.4 GeV, a proposal to 
accommodate it as a 4s state has been made in [6j. Several other suggestions 
have been made for the interpretation of this state, for example, the authors 
of [7j propose it to be a tetra-quark state, others propose a hadronic molecule 
of DiD, DqD* [SI [9], Xci^ [in], XciP [H] and yet another idea is that it is a 
hybrid charmonium [12] or charm baryonium [13], etc. Within the available 
experimental information, none of these suggestions can be completely ruled 
out and its not clear if y(4260) possesses any of these structures dominantly 
or is a mixture of all of them. In Refs. [HI [151 [IS] the authors call the atten- 
tion of the readers to the presence of the opening of the D^D* channel very 
close to the peak position of the y(4260) in the updated data from BABAR 
[TT] and associate the peak corresponding to F(4260) to a D*D* cusp. A fit 
to the data from [U [T7| has been made in [15] and additional presence of a 
rather broad bump around 4.35 GeV has been proposed. 

There are some peculiarities in the experimental findings which motivate 
us to carry out a study of the J/ipiTTT system. There is no enhancement 
found around 4.26 GeV in the process with the D*D* [T^ or other hadron 
final states [I9] and it is concluded that F(4260) has an unusually strong 
coupling to the TiirJ/ip final state [H [2], [3l H]. Further, the data on the 
invariant mass of the vrvr subsystem obtained by the Belle collaboration [Ij, for 
total energy range, 3.8-4.2 GeV, 4.2-4.4 GeV and 4.4-4.6 GeV, have curious 
features. The vrvr mass distribution data in 3.8-4.2 GeV and 4.4-4.6 GeV 
seem to follow the phase space, however, that corresponding to the 4.2- 
4.4 GeV total energy differs significantly from the phase space and shows 
an enhancement near m,r7r = 1 GeV. Do these findings indicate that the 
F(4260) has a strong coupling to fo{980)J/ip, similar to the X(2175) to the 
0/o(98O) ^ 121]? It is interesting to recall that the X(2175) was found 
as a dynamically generated resonance in the ipKK system [221 [23] with the 
KK subsystem possessing the characteristics of the /o(980). Similarly, the 
F(4660) [21] has been suggested as a ■?/''/o(980) resonance [25]. In order to 
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find an answer to this question, we have solved the Faddeev equations for 
the J/ipTiTi and J/ipKK coupled channels and we discuss the formahsm and 
results of our study in the following sections. 

2 Formalism 

In our earlier study of the (pKK system we found the dynamical generation 
of the X(2175) resonance [22]. The study was carried out by solving the 
Faddeev equations for the three-meson system using chiral Lagrangians for 
interaction of the constituent mesons. There are some similarities between 
the X(2175) and the y(4260). Both resonances are of J^^ = 1 nature. 
The X(2175) was found in the 0/o(98O) cross sections [201 [21] and a study of 
this system using chiral dynamics required calculations for the (pKK system 
since the /o(980) is basically a KK molecule in such a formalism. The 
y (4260) has been found in a system of a vector and two pseudoscalar mesons, 
with the two pseudoscalars interacting in the S-wave and with their invariant 
mass showing a dominant peak around 1 GeV in the y(4260) region. This 
hints towards a possibility of clustering of the two pions to form the /o(980). 
If the two pions rearranged themselves to form the /o(980) resonance, the 
F(4260) would be about 200 MeV above the J/'ipfo threshold just as in case of 
the X(2175) which is about 200 MeV above the 0/o(98O) threshold. Besides, 
the diagonal term of the potential obtained from chiral Lagrangian for J/ ijj7^ 
is zero just as the one for the (j)K (or (pir) interaction. However, the 07r (or 
(f)K) scattering matrix is non-zero due to loops of the non-diagonal (coupled 
channel) (pTi — t- KK*{KK*) terms. Similarly, the J/ipir — )■ J/ipTi at the 
lowest order is null but the scattering matrix is formed through iterations of 
the potential involving non diagonal transitions within the coupled channels, 
like J/ipT: — DD* — )■ J/i^tt. This would give rise to three-body diagrams of 
the kind shown in Figs. [1] and [21 All these mentioned similarities between 
the X(2175) and the y(4260), and the experimental findings of y (4260) with 
seemingly strong coupling to the J/ip7nr channel motivate us to carry out a 
three-body calculation of the J/ipirir system. 

We study the J/iprni and J/ipKK systems as coupled channels in the 
isospin base and by considering all the interactions in S-wave. For this we 
solve the Faddeev equations in the formalism developed and used earlier to 
study three-meson system and two meson one baryon systems [221 1211 123 [23 
[221 [30] ■ The different two-meson one baryon systems which we have studied 
so far are the nKN, ttttN and the irKN system, including the corresponding 
coupled channels in each case. In case of S = —1, i.e., ttKN system, we 
find dynamical generation of four isospin 1 and two isospin resonances 
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Figure 1: A three-body interaction diagram where the J /ipir interaction pro- 
ceeds through DD* coupled channel 
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Figure 2: Another possible contribution of the J/ipir amplitude, through 
loops of other coupled channels, to the three-body interaction. 



which we relate to the S(1560), S(1620),S(1660),S(1770), A(1600) and 
A(1810) of the particle data group {PDG) [5j. In 5* = case, we find evidence 
of strong coupling of A^*(1710), iV*(2100), A(1910) to two meson one baryon 
channels [281 129] and also find another isospin 1/2 baryon resonance with 
spin-parity = 1/2"*" around 1920 MeV [2S]. The latter one was earlier found 
in a non-relativistic study of the KKN system using a variational method [21] 
and a signature of this resonance actually seems to be present in experimental 
data on 7j9 — )■ K^K (see Ref. |i32j). We also checked the possibility of 
generation of S* = 1 resonances from three-body dynamics and found a broad 
bump around 1700 MeV [50]. 

The formalism consists of calculating the equation 
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to three-body diagrams with an interaction f which represents a two body 
i-matrix, f = + v^g'^'f, between the jth and kth particles, with j ^ k ^ 
i, and with the ith particle being a spectator. Thus, the superindex on 
the t*'s indicates the non-interacting particle. Such diagrams correspond to 
disconnected three-body diagrams, and removing them one is left with all the 
connected diagrams of the Faddeev equations, the sum of which we denote 
as Tr (where the subscript R denotes the "rest" of the diagrams). Hence the 
matrices correspond to the sum of all the connected diagrams with the 
last two interactions described by f and P . 

These t'-matrices in our formalism are obtained by solving coupled chan- 
nel Bethe-Salpeter equations with the potentials obtained from chiral La- 
grangians. These potentials can be written as a sum of two terms, one de- 
pending only on the center of mass energy of the interacting pair and other 
depending on off-shell variables. As has been discussed elaborately in our 
previous works [22], [26|, [271 EH [291 [30] , interestingly, we find the contribution 
of these off-shell parts of the potentials, and hence of the t-matrices, to cancel 
exactly the three-body forces originating from the same chiral Lagrangian in 
the SU(3) limit. In a realistic case, we found that the contributions from 
the off-shell part of the t'-matrices together with these contact three-body 
terms were negligibly small. An explicit analytic proof of these cancellations 
has been shown for the case of two meson-one baryon systems in [28] and 
for the case of three mesons in [22]. This important finding allows us to 
solve the Faddeev equations by using the on-shell part of the two-body t'- 
matrices, which are independent of the unphysical off-shell parts. Thus, we 
calculate Eqs. ([2|) with the t'-matrices depending on the invariant mass of 
the interacting pair. 



5 



Eqs. ([2]) are solved as a function of the total energy, ^/i and the invariant 
mass of the 23 system, y^S23- In our case, labeling J/ip as the particle 1 and 
KK (and tttt) as the particles 2 and 3, a/s23 is the invariant mass of the two 
pseudoscalars. We define the momenta and other invariant masses in terms 
of a/s and a/s23 as shown in detail in [28] . 

In Eqs. (|2]), the first term, fg^H^ (with i ^ j) represents the simplest 
possible three-body connected diagram which contains two t-matrices, where 
the g^^ matrix elements are three-body Green's functions, defined as 



^2Ek{K + kj)) - E,{k[) - Ej{kj) - EkiK + kj) + it 



(4) 



with Er = \J kr + ml and is the mass of the rth particle of that coupled 
channel to which the element of the matrix corresponds. The mathematical 
expression for the next order diagram, that is the one with three t-matrices, 
is written as fG^^H^g^H'^, where the G^^'^ matrix is 



Qijk _ j 



d^k" 



l9 



ij _ pijk 



(27r)3 

with the elements of the g'^^ being 

1 1 



(*7^J,jV^ = 1,2,3) 



~9''{k",si^) 



2Ei{k") 2Em{k") ^ - Ei{k") - E^{k") + le 
(t^l^m^t^j = 1,2,3) 



(5) 



(6) 



and the matrix F^^'^, with explicit variable dependence, is written as 



r^\k", k', h, s'Z) = t\s09'\k", h) y\k', h) 



-1 r 



tUSr 



(7) 



{j^ry^u = 1,2,3) 

In Eq. (jSD, ^/Slm is the invariant mass of the Im pair which, as mentioned 
above, can be calculated in terms of the variables of the formalism. The 
upper index k" in the invariant mass of Eq. ([7]) indicates a dependence 
with the loop variable as it was shown in [28]. The definition of G^^^ is such 
that the contribution from diagrams with three t-matrices can be written as 

fisiJG'^H^iSru)g'Hkj, h)tHspg) (8) 

which upon substitution of the G*-^ '^ becomes 

= t'(s;j(y ^^^^i^^^'jt^(3.J^7^'(^-,4)t'(5p,) {kj^pj^q = 1,2,3) 
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tHSr 



d^k" ^ ^ r ^ ^ n-1 

T^i9''t^{s';:)g^\k'\k,)[g^\k'^,h)\ 



j^ig''t^{s':i)9'\k\k,)]t\s„,. 



(9) 



(10) 



This shows that in the loops there is a t-matrix with an s^'^ argument re- 
lated to the loop variable k" and the k" integral of all fc" dependent func- 
tions is denoted as G^^^. The next higher order contribution is written as 
^iQijkp Qjki^k gki^i ^ which was numerically shown to be close to that of an 
exact calculation in [28j. In this way Eqs. ([2]) are algebraic equations. The 
remaining inputs required to solve Eqs. ([2]) are the J/V'vr, J /ipK, J/ipK, nn 
and KK t-matrices, the details of which we discuss in the following subsec- 
tions. 



2.1 The t-matrix for the pseudoscalar- vector meson in- 
teraction. 

For constructing the pseudoscalar-vector interaction Lagrangian we follow 
the works in [551 [31] • The starting point for the construction of the La- 
grangian are fields containing all pseudoscalar and vector mesons from a 
15-plet of SU{4:) plus a singlet. In the physical basis these fields read: 
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These two fields differ from tfiose used in [221 El] because of the inclusion of 
a SU{4) singlet in order to take into account the rj-rj' and u-(j) mixing, which 
was not considered in these previous works. 

For each one of these fields a current is defined: 

= {d^^)^-^d^^ (13) 
Jf. = (9^V.)V^ - V.S^V^ (14) 

The Lagrangian is constructed by coupling these currents: 

Cppvv = -^Tr{J^J^). (15) 

The Lagrangian in Eq. (fT5|) is SU{A) symmetric by construction. We 
know, though, that SU{A) symmetry is badly broken in nature, because of 
the heavy charmed quark mass. The first step to break the SU (4) symmetry 
in the Lagrangian is to recognize that the interaction behind the coupling 
of the two currents in Eq. (llSp is the exchange of a vector meson, which 
can be formally visualized within the hidden gauge approach of [351 IMl EZl 
[38] . In this way we suppress terms in the Lagrangian where a heavy meson 
is exchanged by the factor 7 = m|/m^ where rriL is the value of a light 
vector- meson mass (800 MeV) and ttih the value of the heavy vector- meson 
mass (2050 MeV). In the interaction of only heavy mesons {D*Ds, DgD*) 
the vector meson exchanged is the J/ip and we suppress it by the factor 
ip = m\/iv?ji^. We also consider different values for the / appearing in the 
coupling of Eq. ( IT5|) . for light mesons we use / = /tt = 93 MeV but for heavy 
ones f = fD = 165 MeV. 

With our phenomenological Lagrangian we can obtain the potential for a 
given process {P{p)V{k))i {P'{p')V'{k'))j: 

Vij{s,t,u) = -^^{s - u)e.e', (16) 

where s and u are the usual Mandelstam variables, fk is the decay constant 
of the pseudoscalar meson k, e's are the vector-meson polarization vectors 
and i, j refer to the initial and final channels in the coupled channel space. 
The matrix of coefficients can be directly calculated from the Lagrangian 
of Eq. fll5p . Eq. fll6p should be projected into s-wave, which is the only 
partial wave that we study. We come back to technical details in the results 
section. 
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We take the following coupled channels for the strangeness S=l case: 
ir*7r, pK, K*r], K*r]f, uK, (f)K, D*D, D*Ds, J/ipK and K*r]^. And the 
coefficient matrix for these channels in isospin J = | is given below 
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For strangeness S=-l the coupled channels considered are K*7r, pK, K*rj, 
K*-}]/, ujK, (j)K, D*D, D*Ds, J/tpK and K*r]c and the coefficients for these 
channels are the same as for their corresponding S=l channels above. 

For strangeness S=0, one can find the coupled channels and the coefficient 
matrix in [39]. 

To obtain the t-matrix we project in s-wave the potentials of Eq. ( TTGj) 
(removing -e ■ e') and plug them into the scattering equation for the coupled 
channels: 



t = v + vg't. (17) 

In this equation g' is a diagonal matrix with each one of its elements given by 
the loop function for each channel in the coupled channel space. For channel 
i with mesons of masses mi and m2, g'a is given by: 
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where p is the three momentum of the two mesons in their center of mass 
frame. The two parameters ^ and a are not independent, we fix /i=1500 
MeV and use for a the same values used in jS^. These values of a are 
obtained from moderate changes from their natural size |10] in order to fit 
the spectrum for most of the known light and charmed axial resonances. 



2.2 The t-matrix for the pseudoscalar-pseudoscalar me- 
son interaction. 

The TTTT, KK diagonal and non-diagonal potential has been obtained from 
the lowest order chiral Lagrangian j41j 
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The on-shell part of the potential obtained from the Lagrangian Eq. (IT^ in 
S-wave, for total isospin of the two pseudoscalars equal to 0, is 

^KK^KK = 

1 / ml\ 

Kr7r-S>7r7r = 7^ [ ^23 ^T" 1 ' (22) 

These potentials are used to solve the Bethe-Salpeter equations for vrvr and 
KK coupled channels using the same subtraction constants as the ones used 
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in IH] . We would like to mention that we have taken care of the symmetriza- 
tion of the tttt states. The dynamical generation of the a and /o(980) scalar 
resonances in these systems was found using the potentials Eq. fl22|) in [H]. 

Thus, we obtain the t-matrices for the scattering of two pseudoscalars and 
of the vector-pseudoscalar mesons which reproduce the experimental data in 
the corresponding cases. With these inputs we solve the Faddeev equations 
Eqs. 02]). We shall now discuss the results of our calculations. 



3 Results and conclusions 

Using the t-matrices explained in the above section as input, we solve Eqs. 
([2]) for the J/ipnn and J/ipKK channels in total isospin 0, varying the total 
energy ^/s between 4 and 5 GeV and the invariant mass of the two pseu- 
doscalars, ^/s23, between 400 to 1100 MeV. As explained above, the J/ipn 
and J/ tpK interaction is null at the lowest order but it is non-zero when the 
loops of the coupled channels are considered in the iteration of the potential 
leading to the t-matrix. A diagram for the lowest order non-zero contribution 
to the J/ipn interaction has been shown in Fig. [H and its contribution is 
written mathematically as 

VTXTr^-K-Kg^^ Vj/^T,^D*D 9'd* D '^D* D^J /ip-K- (23) 

The potential in Eq. (|T6|) has been obtained by assuming that the momentum 
transfer, i.e., the Mandelstam variable t, in J/ipTi — t- D*D amplitude is 
negligibly small compared to the vector mass. However for the energies and 
channels considered here, such an approximation is not good and we need to 
take the effect of large momentum transfer into account. In order to do this, 
we consider the D*-exchange in the J/ipTi — t- D*D potential (following [12]) 
to get 

r dk' —rrijjt 

Vj/^.^D*D ^ J -^^J/^^^D*DJj^r^^^^^, (24) 

where k' and k are the four vectors of the D* and the J/ip, respectively. This 
would mean that the J/V'tt J/ipn amplitude imphcit in Eq. f l25]) would 
be as shown in Figj3l Similarly, we take into account this correction for the 
J/ipK and the J/ipK amplitudes also. 

With this new potentials we calculate the t- matrix for the J/^-pseudoscalar 
interaction and carry out the calculations for the the J /ipTnr and the J /ipKK 
systems. We find a resonance in both the systems at y/s = 4150 MeV with a 
full width at half maximum of 90 MeV. The peak appears when the invariant 
mass of two pseudoscalars is around that of the /o(980), indicating that the 
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Figure 3: The J /ipir — )■ J /ipir amplitude proceeding through the intermediate 
D*D channel in the loop with a D* exchange at the J/ipn — )■ D*D vertex. 



resonance has a strong coupling to the J/'?/'/o(980) channel. Both the J /ipT^Tt 
and the J/ipKK amplitudes are similar in this energy region, with a differ- 
ence in their magnitudes. We find the J / ipKK amplitude to be much larger 
in magnitude as compared to that of the J/ipTnt system. This reveals the 
strong coupling of the three-body resonance to J/?/'/o(980), since the /o(980) 
couples most strongly to KK [121 El SS] • 

In Fig. m we show the J / tpKK squared amplitude as a function of the 
total energy of the three body system and the invariant mass of the KK 
system. We have also studied the invariant mass spectrum of the two pions 
at \/i = 4 GeV, 4.3 GeV and 4.5 GeV, i.e., in the energy region of the 
resonance and below and above it. To do that we take the three-body 
matrix and multiply it by the phase space factor 

where p = ' — is the momentum of the J/ip in the global center 

of mass system and q = " the momentum of the pion in the 



corresponding two-body center of mass system {Mmv is the invariant mass 
of the two pions) . 

As shown in Figs. O the invariant mass spectrum at ^/s = 4 GeV shows 
a phase space like behavior and the one at y/s = 4.3 GeV shows a dominant 
peak of the /o(980) resonance. At 4.5 GeV, we still see the presence of the 
/o(980) in the two pion mass spectrum but the magnitude of this peak is 
much smaller as compared to the one seen at ^/s = 4.3 GeV, and we find 
that it gradually fades away at higher energies. 

The features described above and depicted in Fig. \5\ agree qualitatively 
with those found for the M^^ spectrum in |1] . One should note that the peak 
of the |Tp matrix is found around 4150 MeV rather than the nominal 4260 
MeV. While 100 MeV difference is not a big difference for a hadronic model 
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Figure 4: [T^P for the J/ipKK system in total isospin / = as a function 
of the total energy, ^/s, and the invariant mass of the KK subsystem, ^/s23■ 



where no parameters have been fitted to the resonance data, the fact remains 
that this difference is the largest one found so far for all the three-body states 
that we have studied before [221 |2S1 ESI ED] • This should be not surprising 
and we would like to attribute it to uncertainties in SU{A) and the fact that, 
unlike other cases, here we have no data to tune our J/tpir and J/tpK (K) 
interaction with our limited freedom in the subtraction constants. 

In order to have some rough estimate of uncertainties we have varied the 
SU(4) symmetry breaking parameter, 7, which enters the evaluation of the 
J/ipTT — )■ J/ipTi or J/ipK{K) — i- J/il)K{K) amplitudes, which proceed as 
shown in Fig. [3] and involve necessarily this parameter. We summarize the 
results here: if we increase 7 in 15 % we find that the strength of the peak 
of Fig. mis also increased in about 50 %. The magnitud of the peaks in Fig. 
EJare also changed in a similar amount. However, we see that the position of 
the peaks and their widths are affected much less and the changes found are 
of the order of 5 MeV for both. 

To summarize the results, the quantum numbers of the state obtained, 
the proximity in the mass to the experimental one and particularly the decay 
mode of the resonance give us strong reasons to associate the state found to 
the y(4260) resonance. 
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Figure 5: |TJjp times the phase space factor for J /ipTm plotted as a function 
of the invariant mass Mtttt of the two pions system for three different total 
energies: a) 4 GeV; b) 4.3 GeV; c) 4.5 GeV. 
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